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ABSTRACT 


Geomagnetic field lines are traced from points on a ballistic 
trajectory to the points where they cross a given altitude shell. 
The intersects thus obtained are called ’’trajectory determined 
intersects.” Three parametric angles identify and describe their 
positions: the look-angles ^and £ , azimuth and elevation respec- 
tively, which relate an intersect to a fixed location on the earth; 
and the line-angle a, which for an observer stationed at this fixed 
location is a measure of the divergence of the line-of-sight from 
the field line. The locus of all trajectory-determined intersects 
is the ’’intersect trace” at the specified altitude level. 

The intersect trace of a sample rocket trajectory is obtained 
for the 100 kilometer level. The variations of the parametric 
angles with time and altitude are calculated for 19 ground loca- 
tions distributed symmetrically around the launch site. For the 
specific flight path considered, there exists an area on the earth 
where optimum parametric angle conditions prevail, i.e., mini- 
mum changes in^ and e, and smallest values of a. 

Finally, the observer-intersect distance is calculated for all 
stations. 
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PARAMETRIC ANGLES OF TRAJECTORY-DETERMINED 
FIELD LINE INTERSECTS 

PREDICTING THE LOCATION OF ARTIFICIAL AURORAS 

by 

E. G. Stassinopoulos 
Goddard Space Flight Center 


INTRODUCTION 

An interesting and promising development in magnetospheric research is the introduction of 
controlled geophysical experiments. 

The first proposal for such an experiment was made by Hess (1965), who suggested that a 
stream of electrons generated by a space-borne accelerator, could produce artificial auroras. 
These auroras would be powerful tools in understanding magnetic field line geometry, conjugate 
point locations, large scale magnetic and electric fields, auroral spectroscopy, plasma instabili- 
ties, and VLF generation. To evaluate the feasibility of the project, Hess suggested in a prelim- 
inary outline of a testing plan (private communication, 1967), that the initial apparatus be flown on 
an Aerobee 350 from Wallops Island. 

In this report we investigate topological aspects of the experiment and appraise the position 
and orientation problems associated With visual or instrumental tracking of artificial auroral 
spots. We have chosen for study a nominal Aerobee 350 trajectory with a launch azimuth of 150 
degrees, launch elevation of 87 degrees, and effective payload of 611 pounds. 

In the context of this study, artificial auroras are auroral streaks generated by a rocket-borne 
accelerator that emits electron beam pulses of constant energy at intervals of three seconds and 
of 1 second duration. The electrons, interacting with the constituents of the upper atmosphere, 
are expected to produce the spots at an altitude which depends largely on their kinetic energy 
(Berger et al., 1968); the more energetic particles penetrating deeper into the atmosphere, the 
less energetic interacting at a greater height, (Figure 17). The sequential appearance of the spots 
in time should follow closely the intersect trace of that altitude level. 

To observe artifieally generated auroras, especially with instruments of relatively small view 
angles, it is desirable to have advance knowledge of the approximate location where the auroral 
spots are expected to appear. If the particle source happens to move along a ballistic trajectory, 
as in the proposed experiment, the observational difficulties are compounded because to nearby 
observers the auroral path may sometimes look erratic. In fact, it will fluctuate strongly in 
azimuth and/or elevation if a station is within a "critical” range around the launch site. Rapid 
displacements of the spots in the field of view of the instruments are then to be anticipated. 
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Sudden and swift variations of the look angles for close-by observers could seriously impair 
the success of a costly experiment. A search -and -aim operation in real time could prove futile 
without some detailed prediction because of (a) the short duration of the. entire experiment, which 
is estimated to last less than five minutes owing to the chosen ballistic flight path, and (b) the 
intermittent very short beam-on intervals. This obstacle can be overcome if the operator antici- 
pates these changes and knows beforehand in what direction they will move. He could, in fact, 
practice the motions in a sort of dry run, long before the real launch. For this purpose a program 
has been designed which can approximately simulate such an experiment. Sample calculations were 
made with a theoretical flight path. Lines of force passing through selected points on this path were 
traced in the direction of the field. Their intersects with a set of successive altitude levels or 
M shells Tf , ten kilometers apart, called trajectory -deter mined intersects ", were determined. 

Three parametric angles identify and describe the position of those intersects: the look- 
angles 0 and e, azimuth and elevation respectively, which relate an intersect to a fixed location on 
the earth; and the line -angle a, which for an observer stationed at this fixed location is a measure 
of the divergence of the line-of-sight from the field line. 

The "intersect trace" or locus of all such intersects for the 100-km shell was obtained. 
Parametric angles for nineteen strategically placed stations were computed. The results were 
plotted versus time and the "critical" range was established. A region of optimum -parametric - 
angle-conditions emerged. 


PARAMETERS OF PROBLEM 

Trajectory, L-band, Time Limitations 

The trajectory used for the computations is within the launch capability of the Wallops Island 
facilities; it pertains to an Aerobee 350 rocket with an effective payload of 611 pounds launched at 
150° azimuth and 87° elevation. Owing to the low latitude of the launch site and the southeasterly 
direction of the impact zone, it crosses relatively few "L-shells". The L-band traversed by the 
flight path is shown in Figure 1; it has a mean value of about 2.56 earth radii. 

Field lines were traced only for points on the trajectory above 150 km. This limited the study 
to the time interval bounded by t T = 110 and t T = 460 seconds, where the prime indicates "after 
liftoff". The evaluation of the parametric angles was further restricted to 200 < t T < 400 seconds 
because during the initial phase of the flight, that is for t T < 200 seconds, unsettled conditions 
prevail. 

Specifically, until its burn-out at about t T = 54 seconds, the rocket is affected by perturbations 
which may cause the actual trajectory to deviate substantially from that originally planned. An 
accurate prediction of the true path to be followed by the rocket cannot be made on the basis of 
initial conditions, such as launch azimuth or elevation, payload, and thrust. Maneuvers and 
activities, required to stabilize the vehicle, orient the apparatus to the field direction, and deploy 
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the electron collecting screen usually would be scheduled during the 146 seconds following burn- 
out. Furthermore, any other task or function that may be necessary prior to the activation of the 
experiment would be performed then. 

During the terminal phase of the flight, when V > 400 seconds, the increasing speed of the 
descending vehicle hampers effective tracking of the spots with instruments of narrow view angles. 
Therefore, with regard to parametric angles, the best time for conducting experiments lies between 
200 and 400 seconds after liftoff. 


Stations 

Figure 2 depicts the location of the nineteen stations selected for investigation; it also shows 
the launch site, the horizontal projection of trajectory points 10 seconds apart, and the correspond- 
ing intersects at 100 km altitude; furthermore, the position of the peak altitude point is indicated 
and the beginning of the "prime" time interval for experiments is marked. 

At first, eleven reference stations were placed in a circle around the peak altitude intersect. 
This arrangement, yielding more information on the look-angle-to-position relationship than any 
other pattern, also saves valuable computer time. Later, eight more stations were added in a 
linear formation to investigate the parametric angles along a line parallel to the main segment of 
the intersect trace. 

METHOD 

Field Line Tracing and Altitude Intersects 

A method developed by the author for geomagnetic field line tracing was presented in a pre- 
vious paper (Roederer et al., 1966). A brief review of this method is given here in Appendix A. 

The procedure used in this report is essentially the same, with only minor modifications to adapt 
it to the present work.* 

The magnetic field strength and field vector components are calculated with McHwain’s new 
MAGNET subroutine (part of the INVAR subroutine), employing the 99-term geomagnetic field 
model by Hendricks and Cain (1966) for the epoch 1960.0, updated to 1965.0. * The L-parameter 
is computed with Hassit and McIlwain T s new (1967) INVAR subroutine. 

The field line integration stepsize, a constant of the tracing process, was set to 2 kilometers; 
tests with larger elements of arc had produced disproportionate errors in the position of the 


*The following changes were introduced into the line tracing procedure: The direction of the tracing was not reversed, since no con- 
jugate intersects were needed; the quantities ' 'arclength’S "dip”, and "declination” were of no concern and were not computed; no 
search for the existence and location of a minimum-B point was made; the L computing part was retained without the value test; the 
section that calculates the intersects was extended to permit their calculation at several consecutive altitude levels during each 
tracing operation. 

tThis model is best known as the GSFC 9/1965. It is available from the Rational ^jpace Sciences Data Center at Goddard Space Flight 
Center, Greenbelt, Maryland 20771. 
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intersects, up to ±.5 degree in latitude and/or longitude, while still expending almost the same 
computer time. 

Lines were traced in the direction of the field only. 


The Look-Azimuth ip 

The azimuth of a point c in space relative to a position a on a sphere is the angle between the 
meridian plane through a and the plane S(a, b, 0), measured clockwise from the northern direction 
of the reference meridian (Figure 3a). The angle 0 in Figure 3b is the azimuth of c with respect 
to a. From spherical trigonometry: 

( cos 6 c - cos 6 cos p\ 

sin 0 sin - p ) = Mp. e c ) • (1) 

For the geoid when a and b are related to the known vectors R and R c (Figure 3a), 0, 9 c 9 and p 
are constants. The cosine of p can be readily found from: 


R * R c 

cos p “ RR ” sin 9 sin 0 c cos [p - 0 c j + cos 6 cos 6 c , (2) 

where 4> and <£ c are the longitudes of R and R c . Thus, p is a function of the four position angles 
of a and b . The azimuth of a field line intersect is easily determined from Equations 2 and 1. 


The Look-Elevation € 

The elevation € of a point in the sky with respect to an observer on the surface of the earth is 
the angle formed by the horizontal plane K and the line from the observer 0 to the point c as in 
Figure 4a. In a geocentric spherical coordinate system it is convenient to evaluate the elevation 
vector ially, in terms of the variables R and R c and the angle p between them. 

Let Q be some vector along the projection of C in the horizontal plane K, pointing away from the 
observer (Figure 4a). Obviously, in unit vector notation, 


sin e = C • R , (3) 

cos e - C * Q . (4) 
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We see from Figure 4b that 


C 



and from Figure 4c and 4d that 


( R x R c ) x R 

IqI 


These relations give: 


1 (r c cosp-r) 

1 (]f[ Rc sin ^) • 


(5) 


(6) 


(7) 


The Line-Angle a 

In the southern magnetic hemisphere, the angle formed by a descending field line B and the 
line of sight c of an observer will be called the "line angle" and will be represented by a (Figure 
5a). Obviously 


„ c • B 

COS p - r^g - “ cos a . 


(8) 


Since 


C = R c -R = f(R Cx -R x )+ ^( R c y - R y ) + k ~( R =,- R 0 


(9) 


the cartesian components of c are 


C x = R c - R x - R c sin 6 Q cos 4> c - R sin 6 cos <f> , 


C = R - R - R sin<9 sin<£ - R sin 6 sin <j> , 

y c y c c r c 


(10) 


C z = R c - R z = R c cos0 c - R cos 9 . 
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Since a unit vector conversion in Equation 9 leads to 


c r = c x sin 6 c cos <f> c + C y sin# c sin<£ c + C z cos(9 c , 
C e - C x cos 0 c cos 4> c + C y cos 6 c sin<£ c - C z sin& c , 

C 0 = - C x sin ^c + C y COS ^c * 

the spherical components are 




^ (ID 


J 


C r 

= c • f 

” R C - R cos p , 


C, 

= C * 6 

- ~ R j sin 6 cos cos(<^ c -<£) - cos 6 cos <9 c 

(12) 


<~e~ 

IU 

il 

- R sin <9 sin [4> c - 



Recalling that 

B r r+ B ^^ + B 0 0 = B = W (l 3 ) 

where the scalar potential v is given by the well known function 

v = 2Z ("f)" + 1 P n m (^) (C cos< ^ + h n m sin^) , ( 14 ) 


we obtain the components of B : 


dv 

d r 


l_dV 

r d6 


1 dW 
r sin 6 d<£> 


(15) 


6 



for increasing r, 0, and<£. The coordinate system used is illustrated in Figure 5b. The field com- 
ponents along the unit vectors u v = - f , u N = - § and u E = are 


B v - - B r : vertically downward 


Bn : horizontally north 


(16) 


B E - B^ : horizontally east 


The field calculations are performed with McHwain’s computer program MAGNET, using the 
99-term Hendricks and Cain model for 1960.0 (GSFC- 9/1965) updated to 1965.0. From the value 
of the field strength and the geometric considerations above, an accurate evaluation of Equation 8 
is possible. 


RESULTS 

The results of our computations are given in Figures 6, 7, 8, 9, 10, and 11 in terms of par- 
ametric angles vs. time. Figures 12 and 13 give the spot- to- station distance with time. For a fixed 
altitude H, the look- elevation e is strictly a function of the line-of-sight distance D: 

e = arc cos (l-H^b 2 ) 1/2 

where D = |c| . Figure 14 is a plot of € vs. D for h = 100 km. It may be used for quick deter- 
mination of the approximate elevation, if the distance is known or vice versa. Predictably, when- 
ever e is greater than 88°, the look-azimuth */> becomes meaningless as a parameter for tracking. 
This condition exists when the intersect trace passes very close to, or over, a station. The effect 
on \p of a near zenith passage of the spots is an abrupt jump in quadrants; it poses problems in 
adjusting instruments that cannot be moved through their local zenith. Stations No. 6 and 7 are 
affected most, whereas Stations No. 1,2,8, and 9, located well to the north or to the south along 
the extended projections of the main intersect segment, experience minimum variations in 0. The 
consequence of these conditions is reflected, for the combined data of e and 0 from all stations, 
in Figures 15 and 16, which depict respectively '’optimum' 1 and "critical" look-angle areas on the 
ground. 

Corresponding to every rocket-generated intersect trace there exists an area on the ground 
where the auroral streaks will be brightest for most of the duration of the flight above 150 kilom- 
eters, the luminance at a given location being proportional to a factor k = 1/sin a . Figure 18 
shows the value of the luminosity factor over the range 0.5° < & <90°. At the observer's magnetic 
zenith the aurora will appear as a faint star, rather than a streak. 
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Using the information available from Figures 10 and 11, the optimum o?-area was determined 
on Figure 19 for k > 2 and k > 3. The apparent brightness of the spots may be of interest to some 
experimenters when selecting the site for a camera station. 

From the evaluation of Figures 15, 16, and 19 the area of "Optimum Parametric Angle Condi- 
tions" emerges and the "Critical Range" develops, as shown in Figure 20. In the "Optimum Area", 
I corresponds to a luminosity factor of k > 2 and IE to k >3. 

Finally, a station- centered polar look-angle plot is added as in Figure 21. It represents the 
intersect trace for each station plotted in e vs. xjj . 
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Appendix 


Field Line Tracing 


The coordinate system used in the field line tracing process is illustrated in Figure 5b. The 
field components along the unit vectors u v , u N , and u E are: 


B : vertical down 


B„, : horizontal north 


B e : horizontal east 


The equations defining an infinite stimal portion of a line of force are: 


SS 

B 


Sr 

B 


rhk _ r cos XS 4> 


B w 




and SS is the element of arc, with components - Sr, rSX. f and r cos \S4> along the unit vectors. The 
coordinates of a generic point of a field line originating at r s , x.,, and <£ s are then given by: 


r r + cr 


\ - k ~ cr 


n 

i = 1 

V 1 

L—j r i B. ss i - 


= 4 >. ~ 


■Dri 

1 E i 

L r i COS ^i B i SSl ’ 


where the sign factor a is ±1 according to whether one proceeds tracing in the same direction as 
the field vector or the opposite. If the origin is in the magnetic southern hemisphere (b v < 0), it 
is necessary to set o < 0 in order to follow the descending line to the nearest intersect. 
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A four point "Adams integration formula” for numerical solutions of differential equations 
is used in the tracing procedure. Its repeated application at equidistant intervals approximates 
the field line. The accuracy of this fit depends on the integration step size. 
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FIGURES AND qAPTIONS 


Figure 

1 Range of magnetic shells crossed by Aerobee trajectory. 

2 Map of stations with trajectory projection and 100-km intersect trace. 

3 Geometry and relationships employed in evaluating the azimuth $ in a spherical 
coordinate system. 

4 Geometry of vectors and angles used in evaluating the elevation e in a spherical 
coordinate system. 

5 Geometry of vectors and angles used in calculating the line angle a. 

6 Variation of look-elevation with time: Stations 1 to 9. 

7 Variation of look-elevation with time: Stations 10 to 19. 

8 Variation of look-azimuth with time: Stations 1 to 9. 

9 Variation of look-azimuth with time: Stations 10 to 19, 

10 Variation of line-angle with time: Stations 1 to 9. 

11 Variation of line-angle with time: Stations 10 to 19. 

12 Variation of distance with time: Stations 1 to 9. 

13 Variation of distance with time: Stations 10 to 19. 

14 Distance dependence of look-elevation for points on the 100-km. altitude shell. 

15 Area of optimum look-angle conditions. 

16 Critical e and \p area. 

17 Approximate interaction altitude of electron beam. 

18 Values of luminosity factor over the range: .5° < a < 90°. 

19 Optimum a area. 

20 Area of optimum parametric angle conditions and critical range. 

21 Station- centered polar look-angle plot. 
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Figure 2— Map of stations with trajectory projection and 100-km intersect trace 









TIME AFTER LIFTOFF (sec) 

Figure 6-Variation of look -elevation with time: Stations 1 to 9. 
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Figure 9— Variation of look-azimuth with time: Stations 10 to 19. 
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Figure 20— Area of optimum parametric angle conditions 
and critical range. 
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